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A new method for separating intensity variations of a source's radio emission hav- 
ing various physical natures is proposed. The method is based on a joint analysis of 
the structure function of the intensity variations and the asymmetry function, which 
is a generalization of the asymmetry coefficient and characterizes the asymmetry of 
the distribution function of the intensity fluctuations on various scales for the inho- 
mogeneities in the diffractive scintillation pattern. Relationships for the asymmetry 
function in the cases of a logarithmic normal distribution of the intensity fluctu- 
ations and a normal distribution of the field fluctuations are derived. Theoretical 
relationships and observational data on interstellar scintillations of pulsars (refrac- 
tive, diffractive, and weak scintillations) are compared. Pulsar scintillations match 
the behavior expected for a normal distribution of the field fluctuations (diffractive 
scintillation) or logarithmic normal distribution of the intensity fluctuations (refrac- 
tive and weak scintillation). Analysis of the asymmetry function is a good test for 
distinguishing scintillations against the background of variations that have different 
origins. 

1. INTRODUCTION 

In analyses of interstellar scintillations of intrinsically variable radio sources, the problem 
of separating variations of the source emission with different origins arises. In particular, 
this is true of separating intrinsic variations of pulsar radio emission, interstellar scintilla- 
tions, and additive and modulational noise. A similar problem appears in the analysis of 
rapid variability of extragalactic sources. In [1-5], it was proposed to use the asymmetry 
coefficient of the distribution function of the radio-emission intensity fluctuations as a test 
for determining the nature of these variations. The asymmetry coefficient is a measure of the 
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deviation of the distribution function of the intensity fluctuations from a normal distribution, 
and is defined as 

7 W - (I)) 2 )\ 312 Ml' 2 ' W 

Here, (I) is the mean intensity of the source's radioe emission, and M 3 and M 2 are the third 
and second central moments of the distribution of the intensity fluctuations. In the case of 
scintillations, 7 is positive, and is related to the scintillation index by the linear relationship 

7 = Am, (2) 

where A is a numerical factor and m is the scintillation index, which is defined as 

_ «/ - am _ Ah 
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If we are dealing with an extragalactic source consisting of a compact scintillating com- 
ponent (core) and an extended non-scintillating component (halo), / corresponds to the 
flux from the scintillating component. The factor A depends on the form of the turbulence 
spectrum and the scintillation regime. In particular, in the case of weak scintillations of a 
pointlike source in the Fraunhofer zone relative to the outer turbulence scale (i.e., the char- 
acteristic size of the largest inhomogeneities), the flux fluctuations are distributed according 
to a Rice-Nakagami law [6, 7], and A = 3/2. In the case of weak scintillations of a point- 
like source in the Fresnel zone relative to the inner scale (i.e., the characteristic size of the 
smallest inhomogeneities), the flux fluctuations are distributed according to a lognormal law 
[8], and A = 3. As is shown in [1], the relationship (2) is a good test for isolating interstellar 
scintillations that give rise to flux variations of extragalactic radio sources. 

However, the asymmetry coefficient is an integrated parameter, which describes the 
summed fluctuations on all scales. Interstellar scintillations can dominate on some time 
scales, while intrinsic variations or noise can dominate on others. To analyze the nature 
of the fluctuations on various scales, it is desirable to introduce some function that would 
yield the value of the asymmetry coefficient for a given time scale in place of the asym- 
metry coefficient. Below, we will introduce such an asymmetry function, 72 (t), which is 
a generalization of the asymmetry coefficient, in the same way that the structure function 
is a generalization of the variance of the temporal process. We will also investigate the 
application of the asymmetry function to the analysis of interstellar scintillations of pulsars. 
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2. DEFINITION OF THE ASYMMETRY FUNCTION 

Let us introdcue the first and second differences 

A 1 (r) = /(* + r)-/(t), (4) 
A 2 (r)=I(t + r)-2I(t) + I(t-r), 

where t is the time and and r is the time shift from instant t. The structure function of the 
flux fluctuations is defined as 

D t (t) = ([^(r)] 2 ) = 2[(/ 2 ) - (I(t + r)I(t))} (5) 

and the mean square of the second difference is 

([A 2 (r)] 2 ) = ({[/(t + r)-/(t)] (6) 
-[I(t)-I(t-r)]} 2 ) = ({[I(t + r)-I(t)] 2 
+ [/(t)-/(t-r)] 2 -2[/(t + r)-/(t)] 
x[I(t)-I(t-T)]})=W I (r)-D I (2T). 

We also introduce the third moment of the second difference A 2 (r) using the formula 

([A 2 (r)] 3 ) = -6(/ 3 ) + 6(/ 2 (t + r)/(t-r)) (7) 
- 12</(t + r)I(t)I(t - r)> + 12(I(t + r)/ 2 (t)> 

and define the second-order asymmetry function as 

72(T) = -([A 2 (r)] 2 )3/- ^ 

The definition of the function 7 2 (r) formally corresponds to the definition of the asymmetry 
coefficient for a random quantity, A 2 (r). Since the third moment of A 2 (r) is negative in 
the case of scintillations, we have added a minus sign in (8), in order for the function 7 2 (r) 
to be positive for intensity fluctuations due to scintillation. 

We also introduce the somewhat different asymmetry function 

72l(r) = _ 2 SMlffl (9) 

72 ' U ) ([A 2 (r)]2)([A 1 (2r)] 2 )V2 ^ 

_ 2 ([A 2 (r)] 3 ) 



[4£) / (r)- J D / (2r)][£> J (2r)]V2- 
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As will be shown below in Appendices A and B, in the case of intensity fluctuations dis- 
tributed according to a Rice-Nakagami law or lognormal law, the function 72,1 (t) for small 
values of r is related to the function [Dj(2t)] 1 ^ 2 / (I) by a formula of the form (2): 

, x A\D i (2t)] 1 / 2 , x 

72,i(r) = [ n {I) n • (10) 

In general, the factor A is a function of r, and we introduce the function A 2; i(r): 

The form of the functions 72 (r), 72,1 (t), and A 21 (r) for the two distribution laws (a normal 
distribution for wave-field fluctuations and a lognormal law for intensity fluctuations) is 
determined in Appendices A and B. 



3. WEAK SCINTILLATIONS 

We will consider the properties of the asymmetry function in the case of scintillations 
of radio sources on electron-density inhomogeneities of plasma with a power-law turbulence 
spectrum. It is known that, in this case, the scintillation regimes are determined by the 
structure function of the phase fluctuations on a scale that is equal to the scale of the first 
Fresnel zone Ds(bF r ), where 6p r = (feff/k) 1 ^ 2 , r e ff is the effective distance from the observer 
or the source to the layer containing the turbulent medium (in the case of a statistically 
uniform medium, this is the distance from the observer to the source), and k = 2tt/\ is the 
wavenumber. If 

D s {b Fr ) < 1, (12) 

the scintillations are weak, the characteristic spatial scale giving rise to the intensity fluctu- 
ations is equal to 6p r , and the scintillation index can be found from the formula [9] 

m 2 = K(n)D s (b Fr ), (13) 

where K{n) is a numerical factor of the order of unity that depends on the exponent n of 
the turbulence spectrum. 

Numerical calculations of the one-dimensional distribution function of the flux fluctua- 
tions F(I) were carried out in [10, 11], and the second and third moments of the intensity 
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Figure 1. Time variations of the emission intensity for PSR B1933+16 at 610 MHz 
(top), PSR B1642-03 at 4.85 GHz (center), and PSR B0329+54 at 4.85 GHz (bottom). 

fluctuations in the case of a Kolmogorov spectrum for the inhomogeneities of the index of 
refraction for various values of the inner turbulence scale were determined. It was shown 
that, in the weak scintillation regime, as well in the strong scintillation regime with values 
of DsipFr) that are not too large, the distribution function /(/) is close to a lognormal 
function. Using the results of these studies, we find A = 2.78 for a Kolmogorov spectrum 
without a turnover and A = 2.86 for a Kolmogorov spectrum with an inner scale I equal to 
the scale bp of the first Fresnel zone {I /bp = 1). Thus, with increase in the inner scale of the 
turbulence spectrum, the factor A increases and approaches three, which corresponds to a 
lognormal distribution for the flux fluctuations. The calculations show that relationship (2) 
is valid up to a value of m that approaches unity from the side of the weak (unsaturated) 
scintillation regime. 

Note that the one-dimensional distribution function /(/) and asymmetry coefficient 7 
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B1933+16 f= 610 MHz 




Figure 2. Structure function, asymmetry function, and coefficient A2,i 
versus the time shift for refractive scintillations of PSR B1933+16. 

describe the flux fluctuations on the main scale of the spatial pattern of the scintillations and, 
accordingly, on the main temporal scale To of the structure function of the flux fluctuations. 
When r > r , the functions 72 (t), 72,1 (t), and A 2 ^(r) will be determined by asymptotic 
relationships that correspond to a lognormal distribution law for the intensity fluctuations 
(see Appendix B): 

3 £>j(2t) 
72 ^ = 2 [AD^t)- D^t)} 1 / 2 ^ ( ' 

72 ,i(r)=3[ J D / (2r)] 1 / 2 /a), (15) 

^2,i(r) = 3. (16) 

Flux fluctuations with small spatial and temporal scales are determined by the diffraction 
of the radio waves on inhomogeneities in the index of refraction with sizes that are much 
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Figure 3. Same as Fig. 2 for PSR B0329+54. 

smaller than the first Fresnel zone. Therefore, the distribution of small-scale flux fluctuations 
will more closely fit a normal distribution law, and the functions 72 (t), 72,1 (t), and A 2 ,i(t) 
will better match the asymptotic relationships (A.8)-(A.10) of Appendix A. 

4. DIFFRACTIVE SCINTILLATIONS 

For large values of the parameter D s (bp r ), the intensity fluctuations consist of two com- 
ponents: small-scale (diffractive) and large-scale (refractive). With increasing D s (bp r ), the 
distribution function of the field of the diffractive component tends to a normal distribution 
[9], the scintillation index tends to unity, and the characteristic spatial scale of the diffraction 
pattern corresponds to the field coherence scale l/k9 scat , where scat is the characteristics 
scattering angle. When r < r , the functions 72 (r), 72,1 (t), and A 2j i(t) will correspond to 
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the asymptotic expressions (Appendix A) 



72 (t) 



3 Di(2t) 

4 [4D / (r)-D / (2r)]V2(/)' 



(17) 



72,i (r) 



3 [^/(2r)] 1 / 2 
2 (/) 



(18) 




72,i (t) 



(19) 



[ J D / (2r)] 1 /2 2 ' 



and, when r ^> r , will be determined as 72(1") = 2/v^6, 72,1 (t) — 2\/2, A 2 ,i(r) = 2. 



The refractive component is determined by weak focusing of the radiation on inhomo- 
geneities with characteristic scales of the order of the scattering disk: r e ff9 scat [9]. Refractive 
scintillations can be treated as weak scintillations on large-scale inhomogeneities, with sizes 
much larger than the scale of the first Fresnel zone. Therefore, we expect a lognormal distri- 
bution law for the flux fluctuations. When r > r , the functions 72(1"), 72,1 ( r ), and A 2; i(t) 
will be determined by the asymptotic relationships (14)-(16). 



We used data on diffractive, refractive, and weak scintillations of several pulsars to ana- 
lyze the applicability of the above relationships for the asymmetry function to actual obser- 
vations of interstellar scintillations of pulsars. 



We used the five-year monitoring data for PSR B1933+16 and B0329+54 (1990-1995) 
at 610 MHz [12] (NRAO, Green Bank) to analyze refractive scintillations. The observations 
were carried out in a 16-MHz band; diffractive scintillations were considerably smoothed 
by averaging the signal over about an hour of observations. Figure 1 (top) shows the time 
variations of the emission intensity for PSR B1933+16. The time scale of these variations 
is 33 days [12]. We used relationships (5), (9), and (11) to calculate the structure function 



5. REFRACTIVE SCINTILLATIONS 



6. COMPARISON WITH OBSERVATIONS 



Refractive and Weak Scintillations 
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Figure 4. Same as Fig. 3 but for weak scintillations. 

of the flux fluctuations D(t), asymmetry function 72,1(7"), and coefficient v4 2j i(r) for all 
the types of scintillations considered below. Before calculating these functions, we have 
normalized the time series to their average values. 

Figures 2 and 3 show D(t), 72,1(7"), and A 2 ,i(r) as functions of the time shift in days for 
PSR B1933+16 and B0329+54, respectively. The structure function reaches a constant level 
equal to 2m 2 , where m is the modulation index of the intensity variations. The time shift at 
which -D(r) decreases by a factor or two gives the time scale of these variations. In 72,1(7") 
and Ac(t), this time scale corresponds to a shift that is a factor of two smaller [see formula 
(10)]. We can see from Figs. 2 and 3 that A 2 ,i(r) is approximately constant for shifts of less 
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Figure 5. Same as Fig. 4 but for PSR B1642-03. 

than 40 days, and this average level is A 2i i ~ 3 for PSR B1933+16, based on formula (5) for 
refractive scintillations, and A 2) \ ~ 2 for PSR B0329+54. The difference of A 2) i from three 
for PSR B0329+54 may be due to the fact that the index of the electron-density fluctuation 
spectrum toward this pulsar differs from the Kolmogorov index [13]: it is n — 3.5. It seems 
reasonable to suppose that the decrease of the spectral index should result in a decrease in 
A 2i i at small time shifts. 

We used observations of the pulsars PSR B0329+54 and B1642— 03 carried out at 
4.85 GHz on the 100-m Effelsberg telescope for our analysis of weak scintillations. The input 
data and observing parameters are given in [13, 14]. Before recording the signal, we averaged 
the pulsed emission over 10 P x for PSR B0329+54 and 39 Pi for PSR B1642-03 (where P 1 
is the pulsar period). The time resolutions were 7.15 and 15 s, respectively. Figure 1 shows 
the time series for these pulsars (central and bottom graphs). We used the fourth channel 
(the observations were carried out in four frequency channels of 60 MHz each), which was 
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Figure 6. Top: average profile of PSR B1133+16 at 111 MHz. Center: pulse- 
to-pulse intensity variations in a single 1.25-kHz channel (for each pulse the inten- 
sity was obtained by integrating within the longitude interval shown with dashes on 
the average profile). Bottom: autocorrelation function of these intensity variations. 

the most free of interference, for PSR B0329+54. The modulation indices for the intensity 
variations for PSR B0329+54 and B1642-03 were 0.45 and 0.64, respectively. 

Figures 4 and 5 show D(r), 72,1 (t), and A 2j i(r) as functions of the time shift in seconds. 
The maximum time shift for which these functions were calculated was one-fourth of the 
total time interval of the observations. The scintillation time scales for PSR B0329+54 and 
B1642— 03 are 2200 and 1170 s, respectively (at the 1/e level of their autocorrelation func- 
tions). On time scales of 1000-2000 s for PSR B0329+54 and 400-700 s for PSR B1642-03, 
the values of A 2j i(r) are 2.8 and 3.3, respectively. These agree well with the value A 2> i = 3 
predicted by the theory for weak pulsations [formula (16)]. The constant level at the shortest 
time shifts in the structure function D{t) is due to uncorrelated noise. 
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We did not correct the calculated functions for noise, because, as was shown by the 
analysis, their third moments are nonzero, and this correction would introduce strong dis- 
tortions to 72,1 (t) and A 2j i(r). The constant level A 2 ,i ~ 2 at time shifts r < 200 s for 
PSR B0329+54 is probably due to diffractive scintillations, whose time scale at 610 MHz is 
t d « 90 s [13]. For PSR B1642-03, t d « 130 s at the same frequency [14]. For this pulsar 
there are strong fluctuations in A 2 ,i(r) at shifts shorter than 200 s, due to the effect of noise; 
however, the mean level of about 1.5 is probably associated with diffractive scintillations on 
small-scale inhomogeneities. 

Diffractive Scintillations 

We used observations of PSR B1133+16 carried out on September 27 and 29, 2001 on the 
Large Phased Array of the Lebedev Physical Institute in Pushchino at 111 MHz to analyze 
the intensity variations of this pulsar using the asymmetry function. Individual pulses of 
the pulsar were recorded in 96 channels of the receiver (the bandwidth of each channel was 
1.25 kHz) with a time resolution of 1.56 ms during 200 s, and were stored on a computer 
disk. We performed the structural analysis after removing the dispersion shift in all frequency 
channels and integrating the spectrum in the selected interval of pulse longitudes. Figure 6 
(top) shows the mean pulse of PSR B1133+16, obtained by accumulating pulses during 200 s 
in a 120-kHz band. The same figure (center) presents the pulse-to-pulse intensity variations 
in one of the channels (1.25 kHz). The intensity was obtained for each pulse in all the 
channels by integrating the signal in the longitude interval of the first component of the 
mean profile, which is denoted in Fig. 6 with dashes. The same figure (bottom) shows the 
autocorrelation function (ACF) calculated for the intensity fluctuations I(t) in this channel. 
The ACF has three time scales: the first is due to pulse-to-pulse variations within one pulsar 
period (Pi = 1.188 s), the secondm which is about 7 s, is due to diffractive scintillations, and 
the third represents slower variations on a time scale of ~20 s. The fact that the 7-s time 
scale is due to diffractive scintillations follows from the fact that the decorrelation of the 
spectra (96 channels 1.25 kHz each) occurs precisely on this time scale (with pulses separated 
by 6 Pi). To reduce the effect of intrinsic rapid variations of the pulsar, we smoothed the 
input I(t) data in all the channels on a three-point interval (3 Pi). When calculating D(r), 
72,1 ( r ), and A 2;1 (r), we averaged the second and third moments of I(t) over all frequency 
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Figure 7. Structure function, asymmetry function, and coefficient A2,i versus time shift for 
PSR B1133+16 at 111 MHz averaged over two observational days (diffractive scintillations). 

channels. Figure 7 shows the obtained functions averaged over two observational days. We 
can see that the different time scales are well separated in the functions 72,1 (t) and y4 2j i(r). 
The value of y4 2 i(r) on scales < 7 s is 1.4, in good agreement with the theoretical value of 
1.5, corresponding to diffractive scintillations. 

7. CONCLUSION 

Our analysis of observed pulsar intensity variations using the asymmetry function 72,1 (t) 
and function A 2 i(r) has shown that these variations match the behavior expected for a nor- 
mal distribution for the field fluctuations (diffractive scintillations) or a lognormal distribu- 
tion for the intensity fluctuations (refractive scintillations, weak scintillations), in accordance 
with the derived theoretical relationships. For all the scintillation regimes, there is a func- 
tional relationship between the asymmetry function and the structure function of the inten- 
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sity fluctuations. This relationship has the simplest form for the asymmetry function 72,1 (t) 

72 i(t) 

at small values of r: 72,1 (t) oc D 2 ,i{'2t) . Accordingly, the function A 2) i(r) = (I) ' 

V-D/(2r) 

is a constant. This indicates that analysis of the function A 2 ,i(r) could provide a good test 
for isolating scintillations against the background of variations that have other origins, and 
that the proposed method can be efficiently used to study and separate out intensity vari- 
ations that have different physical natures. In particular, it can be applied to the analysis 
of flux variations of extragalactic sources, in order to distinguish variations due to intrinsic 
variability of the source and effects associated with the propagation of the radiation in the 
interstellar plasma. 

We note also that there are prospects for further studies in this direction. It is necessary 
to study the dependence of the form of the asymmetry function on the spectrum of the 
inhomogeneities in the index of refraction, in particular, for a power-law spectrum, and the 
dependence of the asymmetry function on the spectral index. Studies of the dependence 
of the form of the asymmetry function on the source structure are also required. If such 
theoretical calculations, confirmed by experimental data, are available, measurements of the 
asymmetry function can yield additional information about both the medium and source. 
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THEORETICAL RELATIONSHIP BETWEEN THE ASYMMETRY FUNCTIONS 7 2 (r), 
72,1 (r) AND THE STRUCTURE FUNCTION D/(r) FOR A NORMAL DISTRIBUTION 

FOR THE FIELD FLUCTUATIONS 

Let us assume that the wave-field fluctuations have a normal distribution and that the 
average value of the field is zero: 

(E(t)) = 0. (A.l) 

This equality implies that second-order moments of the form (E 2 (ti)) are also equal to zero 
[9]: 

(E\ tl )) = ([E*( tl )}>) = 0. (A.2) 

Formulas (A.l) and (A.2) correspond to diffractive scintillations in the saturated scintil- 
lation regime [9]. In this case, the third moment of the intensity or the sixth moment of the 
field can be represented as a sum of products of second-order field-coherence functions [15]: 

(I( tl )I(t 2 )I(t 3 )) (A.3) 

= (E(t 1 )E*(t 1 ))(E(t 2 )E*(t 2 ))(E(t 3 )E*(t 3 )) 
+ (E(t 1 )E*(t 2 ))(E(t 2 )E*(t 1 ))(E(t 3 )E*(t 3 )) 
+ (E(t 1 )E*(t 3 ))(E(t 2 )E*(t 2 ))(E(t 3 )E*(t 1 )) 
+ (E(t 1 )E*(t 1 ))(E(t 2 )E*(t 3 ))(E(t 3 )E*(t 2 )) 
+ (E(t 1 )E*(t 2 ))(E(t 2 )E*(t 3 ))(E(t 3 )E*(t 2 )) 
+ (E(t 1 )E*(t 3 ))(E(t 2 )E*(t 1 ))(E(t 3 )E*(t 2 )) 
= (I) 3 + (I)[B I (t 1 -t 2 ) + B I (t 1 -t 3 ) 

+ £ 7 (t 2 - t 3 )\ + B E (t! - t 2 )B E (t 2 - t 3 ) 

x B E (t 3 - t 2 ) + Bsih - t 3 )B E (t 2 - h) 
x B E (t 3 -t 2 ). 



Here, 



B E {t l -t k ) = {E{t l )E*{t k )) (A.4) 
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is the field-coherence function and Bi(ti — tk) = ([I(U) ~ — (H^k))]) = Be(U — 

tk) x B* E {ti — tk) is the correlation function of the intensity fluctuations. Using formula 
(A. 3) in the calculation in (7), we obtain a relationship for the third moment of the second 
difference: 



([A 2 (r)] 3 ) = -12(/) 3 + 12(/)5 / (2r) 
+ 24<J)£ / (r) - 2AB e (2t)B i (t) = -6(1) D^r) 



(A.5) 



+ 24(7) (7) 



{I? ~ 2^(2r) 



1/2 



127^(2^(7) 



(7) 2 - -D,(2r) 



1/2 



Let r be the characteristic scale of the structure function of the intensity fluctuations 
Di(t). For small values, r < r , we obtain the asymptotic formula 



<[A 2 (r)] 3 ) = -M2T) 



4D 7 (r) - Dj(2r) 



;a.6) 



4 — (7) 

For large values, r ^> r , using equality Dj{r) = 2((7 — (7)) 2 ) = 2(7) 2 yields 

<[A 2 (r)] 3 ) = -12<7) 3 . (A.7) 

Using these asymptotic relationships for ([A 2 (t)] 3 ) and formulas (6) and (8), we derive 
for the function 7 2 (r) the asymptotic expressions 

3 Dj(2t) 

.MJ./2//A ^ - T °' ^ 



72 (t) 



4 [4£) J (r)- J D J (2r)]V2(/)' 

7!' r>>ro - 



Similarly, we have 



[ J D / (2r)] 1 /V(/) ) T <r , 



72,i (r) 



(A.9) 



2>/2, r > ro . 

We see that, for small values of r, 72,1(1") is related to [D i (2t)] 1 ^ 2 / (I) by the same formula 
as that relating the asymmetry coefficient to the scintillation index. 
We obtain for A 2 ,i(t) 



A 2 ,i(r) = (I) 



72,i (t) 



[^/(2t)]V2 2 
A 2:1 (t) = 2, r>r . 



(A.10) 
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THEORETICAL RELATIONSHIP BETWEEN THE ASYMMETRY FUNCTIONS 7 2 (r), 
72,i (r) AND THE STRUCTURE FUNCTION D 7 (r) FOR A LOGNORMAL 
DISTRIBUTION FOR THE INTENSITY FLUCTUATIONS 

Let us represent the intensity I(t) as 

I(t) = I exp[ X (t)] (B.l) 

and assume that the fluctuation distribution law x(t) is normal. For a normally distributed 
quantity, the following equality is valid [8]: 

1 



Jo(exp[x(t)]) = 7 exp 



(x(t)} + -((s x (t)) 2 } 



(B.2) 



Taking into account the fact that the mean intensity does not change during the propagation 
of the wave in a turbulent medium, i.e., (I) — Io, we obtain 

<*(*)> = ~(mt)) 2 ). (b.3) 

We find with (B.2) and (B.3) 

(J 2 ) = (/ 2 exp[2 X (t)]) (B.4) 
= / 2 exp [2( X (t))+2((5 X (t)f}] 
= / 2 exp [((5 X (t)) 2 )} , 



(I(t)I(t + r)> = exp [(5x(t)5 X (t + r)>] (B.5) 
= (/ 2 )exp 



where 



D x (r) = ([ X (t + r) - X (t)] 2 ) (B.6) 
= 2([5x(t)] 2 )-2(5x(t)5 X (t + r)) 

is the structure function of the fluctuations, x(t)- 

Formulas (B.4) and (B.5) enable us to express the quantities exp [((#xW) 2 )] an d 
exp [— |D x (r)] in terms of the scintillation index and the structure function of the intensity 
fluctuations: 

exp[((5x(t)) 2 )] = l + m 2 , (B.7) 



Then, using the formula 



(iitMtjifo)) 

J 3 exp \3{(S X ) 2 ) - \D x {h - t 2 



/ 3 (l + m 2 ) 



x exp 



1, 



--[D^-^ + D^-h) 



+ D x {t 2 -h)\ 



we obtain for the third moment of the second difference ^(t) 

<[A 2 (r)] 3 > = -6/ 3 exp [3((5 X ) 2 )] 
x |l - exp[-D x (2r)] -2exp[-D x (r) 

+ 2 exp 

Using (B.4) and (B.5), we find 



-D x (t) - \d x {2t) 



2\3 



<[A 2 (r)] 3 > = -6/ 3 (l + m 2 ) 



x 1 



x 



2/ 2 (l + m 2 ). 



+ 2 



/ 2 (l+m 2 ) 



2/ 2 (l + m 2 ) 

For small m 2 < 1, we obtain the asymptotic relationship 



<[A 2 (r)] 3 > = -D I (2r) 



ADj{t) - Dj(2t) 



2 I 
Substituting this asymptotic relationship into (8)-(10) yields for m 2 



72 (r) = - 



D t (2t) 



2 [4 J D / (r)- J D/(2r)]V2/ ' 



72 , 1 (r) = 3[ J D / (2r)] 1/2 //o, 
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